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In this work, I report the latest lattice QCD calculations of nucleon and hyperon 
structure from chiral fermions in 2+1-flavor dynamical simulations. All calculations 
are done with a chirally symmetric fermion action, domain-wall fermions, for valence 
quarks. I begin with the latest lattice results on the nucleon structure, focusing 
on results from RBC/UKQCD using 2+1-flavor chiral fermion actions. We find the 
chiral-extrapolated axial coupling constant at physical pion mass point to be 1.23(5), 
consistent with experimental value. The renormalization constants for the structure 
Q ' functions are obtained from RI/MOM-scheme non-perturbative renormalization. We 

^ : 

find first moments of the polarized and unpolarized nucleon structure functions at 

^ , zero transfer momentum to be 0.133(13) and 0.203(23) respectively, using contin- 

I uum chiral extrapolation. These are consistent with the experimental values, unlike 

OO ■ 

■ previous calculations which have been 50% larger. We also have a prediction for the 



transversity, which we find to be 0.56(4). The twist-3 matrix element is consistent 



o 

^ \ with zero which agrees with the prediction of the Wandzura-Wilczek relation. 



In the second half of this work, I report an indirect dynamical estimation of the 
strangeness proton magnetic moments using mixed actions. With the analysis of 
hyperon form factors and using charge symmetry, the strangeness of proton is found 
to be —0.066(26), consistent with the Adelaide- JLab Collaboration's result. The 
hyperon S and H axial coupling constants are also performed for the first time in a 
lattice calculation, = 0.441(14) and ^hs = -0.277(11). 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) has been successful in describing many properties of 
the strong interaction. In the weak-couphng regime, we can rely on perturbation theory to 
work out the path integral which describes physical observables of interest. However, for 
long distances perturbative QCD no longer converges. Lattice QCD allows us to calculate 
these quantities from first principles. 

In Lattice QCD, space and time are discretized in a finite volume, and we use Monte 
Carlo integration to evaluate the remaining integral. Since the real world is continuous and 
infinitely large, at the end of the day we will have to take a —>■ and V oo. However, 
using current computer resources, we cannot yet simulate full QCD at the physical pion mass. 
With the help of the chiral perturbation theory and calculations at multiple heavier pion 
masses which are affordable in terms of available computational resources, we can extrapolate 
quantities of interest to the physical limit. Such calculations also help to determine the low- 
energy constants of the chiral effective theory. 

There are a few choices of fermion action that have been commonly used in lattice QCD 
calculations. Each has its own pros and cons. They differ primarily by how they maintain 
symmetry, their calculation cost and their discretization error. The consistency of results 
from different fermion actions demonstrates the university of discretizations from lattice 
QCD. 

The most expeusive Cass of disc.et.ati„n a.e the ch.al fe™io„ aetionsfl fl 

domain-wall fermion (DWF) or overlap fermion. Such actions maintain the chiral symmetry 
of the fermions at great cost, but for this cost, we derive significant benefits. The theory 
is automatically 0{a) improved which makes it particularly well suited for spin physics and 
weak matrix elements. Since symmetry remains at the non-zero lattice spacing, it further 
simplifies the renormalization calculation on the lattice and the chiral extrapolation. 

A much cheaper alternative is the (improved) staggered fermion action (asqtad) 
The relatively fast simulation has great potential to be the first lattice calculation to reach 
the physical pion mass with 2+1 flavors. However, it introduces the problem of taste: 
each single fermion in the action contains four tastes. Although these extra tastes can 
be removed by the "root trick", mixing among different parities and tastes remains in the 
theory. There have been spirited debates over whether it is proper to use asqtad j^, 9|, but for 
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practical purposes no evidence of anomalous results have yet been found; theoretical proof 
of its correctness is an ongoing effort. However, the issue of taste-breaking makes baryonic 
operators a nightmare to work with, regardless of any potential cost savings. 

One might combine the best features of both of these actions by using a mixed action. 
Cheap staggered fermions can be used for the expensive sea quarks, while chiral domain-wall 
fermions are used in the valence sector, where they protect operators from mixing. We will 
discuss the employment of such an action in the final section. 

In this work, we use lattice QCD techniques to pursue long-distance physics, solving 
non-perturbative QCD from first principles. The structure of this article is as follows: 
In Sec. ini I report the latest 2+lf DWF valence calculation using RBC/UKQCD gauge 
configurations. In this work, we concentrate on the results for the axial coupling constant, 
first moment of the unpolarized quark and helicity distributions, and the zeroth moment of 
the transversity and twist-3 matrix element di. We reproduce the experimental numbers 
for the first three quantities using continuum perturbation theory, predict the value of the 
transversity at leading order and check the Wandzura-Wilczek relation. We also compare 
our results with those of other lattice groups. In the Sec. IIIIl I present some work done by 
the Jefferson Lab hyperon project, calculating the proton strange magnetic moment and the 
axial couplings of the hyperons. This is the first time that the strange magnetic moment is 
done using dynamical lattice data, and we find a value consistent with the Adelaide- JLab 
Collaboration's result. The S and S axial coupling is for the first time done in lattice 
QCD. We find the numbers are more accurate than what had been estimated by the chiral 
perturbation theory or large- iVc theory. 



II. NUCLEON STRUCTURE 
A. Lattice Parameters 



In this calculation, we uses lattices generated by the RBC/UKQCD collaborations with 
the chiral DWF action and a full 2-|-l-flavor dynamical quark sector. The ensembles range 
in pion mass from 625 down to 300 MeV, at a lattice scale 1.6 GeV in a 3 fm box. The 
details of the gauge configurations can be found in Ref. 10[ for hadron properties in a 2 fm 
box. 
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TABLE I: 2+1 flavor Gaussian-smeared source parameters 



'^sca 


0.005 


0.01 


0.02 


0.03 


^src 


0,32 


0,16,32,48 


0,16,32,48 


0,16,32,48 


^ of conf. 


180 


119 


49 


54 


(GeV) 


0.319(3) 


0.399(3) 


0.535(3) 


0.625(3) 


niN (GeV) 


1.085(16) 


1.169(19) 


1.204(13) 


1.474(18) 



On these ensembles, we use a Gaussian-smeared source to improve the signals. The 
source-sink separation is fixed at 12 time units. The choices of sink and the corresponding 
number of configurations in this work are listed in Table [B 

The interpolating field used in our calculation is 

X"" =Y. [ul{y,,t)C^,db{y2,t)\ ne,«(y3,t)0(yi - x)0(y2 - x)(t>{y, - x) (1) 

X ,a,b,c 

We calculate the nucleon two-point function (C2pt of with smearing parameters A and 
B) as 



(2) 



and three-point function which is defined as 



AB 



(rl ^('^sro ti '^snk)) 



src) P src)), 



(3) 



where O is 
Refs. [n 



;he operator of interest. For more details, please refer to our earlier work in 
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13|. 



When one calculates the three-point Green function, there are two possible contraction 

igure[Tl Disconnected 



17|. It would 
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topologies: "connected" and "disconnected" diagrams, as depicted in 
pieces are notoriously difficult to calculate directly on the lattice [ 
require numerous source vectors in the fermion-matrix inversion, and there was no reasonable 
way of calculating these contributions when this calculations started. However, this difficulty 
might be resolved in the near future with development of new techniques by members of the 



FIG. 1: The three-point Green function contains connected (left) and disconnected (right) dia- 
grams. The disconnected piece is ignored in this work by focusing only on isovector quantities. 



USQCD collaboration. In this work, only "isovector" quantities will be discussed, in which 
the disconnected piece cancels under isospin symmetry. 



We apply nonperturbative renormalization (NPR) in RI/MOM scheme ISf] to the above 
quantities. In general, the operators of interest can mix with lower-dimension operators as 



Oi{n) = Zi{n, a)Oi(/i) + ^ %(/i; a)Oj{n) 



(4) 



With the good chiral symmetry of DWF action, however, we are protected from this mixing 
problem. We calculate Zo^if^] O') in RI/MOM scheme, where /i must fall inside the renormal- 
ization window Aqcd ^ ^ 1/a. Then we convert to MS scheme 19|, running to 2 GeV 
to get the renormalization constant for the operators. In this work, the NPR is done on a 
smaller lattice ensemble, the 2 fm ones; since it is a short-distance quantity, the renormal- 
ization constants are not as sensitive to finite-volume effects as other observables. Detailed 
step-by-step descriptions can be found in Refs. 
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121 ]. The renormalization constant in 



RI/MOM and MS scheme for first moment of the unpolarized distribution and helicity and 
zeroth moment of the transversity (from top to bottom) are shown in the Figure [2l 



B. Numerical Results 

1. Axial Coupling Constant 

The axial charge is well measured in the neutron /9 decay experiment and hence it is 
a natural candidate to demonstrate how well the lattice QCD approach with the chiral 
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FIG. 2: Renormalization constants in the chiral limit. The hghtly- filled circles are the renormaliza- 
tion constants in RI/MOM scheme, and the diamond points are MS-scheme at 2 GeV. The fitted 
lines are used to remove (ap)^ artifacts. 

extrapolation to the physical pion point. The isovector vector and axial charges, gv and qa, 
are defined as the zero-momentum-transfer limits of the following. 



{p\V^^{x■, q = 0)|n) = Up {-fi^gv) Un 

= 0)|n) = Mp (7m75^a) Wn- 



(5) 
(6) 



Because of chiral symmetry on our fermion action, a Takahashi-Ward identity ensures that 
the two currents, which are related by chiral transformation, share a common renormaliza- 
tion: Za = Zy up to a lattice discretization error of 0{ma?'). Since the vector current is 
conserved, its renormalization is easily obtained as the inverse of the vector charge gy- Thus, 
by calculating the ratio of the three-point functions for gA/ gvi we get the renormalized axial 
charge, [gA)^'^^- The results are shown as red triangles in Figure [31 

In order to reach the physical pion mass value for our result, we adopt the chiral ex- 
trapolation expression from the the small-scale expansion (SSE) scheme 20|]. In this scheme, 
one uses explicit degrees of freedom from the pion, nucleon and A(1232), then expands in 
terms of Aq, the mass splitting between the and A in the chiral limit, which is treated 
as 0(e). Here we adopt a formulation which is correct up to O(e^), as seen in Ref. 2l|, |22| . 
and we also try its finite- volume corrected form. The grey band in Figure [3] indicates the 
uncertainty due to the jackknife chiral extrapolation with SSE formulation; our preliminary 
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extrapolated axial coupling constant is 1.23(5), consistent with experimental data. If we 
add finite-volume effects, with the lattice box fixed at 3 fm, we can see the correction at 



our pion mass points are tiny, as indicated by the green band. In_^right of Figure [3], we 
compare our result to that of other lattice groups 2J, uA 



and to our previous 



calculation Il3l]. Compared with our previous DWF quenched 2.4 fm and 2-fiavor 2 fm 
box calculations, we see a clear consistency in the axial charge coupling. In fact, most of the 
calculations (either using chiral fermion or not) are consistent with each other. One thing 
to notice is that we currently have the lightest dynamical pion mass point, although the 
result from the lightest ensemble (in both our 2- and 2+1-flavor) deviates from the chiral 
extrapolation curve. We have more gauge ensembles in the 2-|-l-fiavor case, and we will be 
able to verify in the near future whether this is just simply due a lack of statistics or to a 
finite-volume effect that is not accounted accurately by the chiral perturbation theory. 
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FIG. 3: (left) Renormalized axial charge in terms of pseudoscalar mass with the small scale expan- 
sion fit (grey band) and an estimation of the finite- volume effect (green band) 
(right) The axial charge comparison from various lattice group 



2. Unpolarized quark and helicity distribution 
The moments of the unpolarized quark distribution and helicity distribution are defined 



as: 



c/xx"[g(x) - (-l)"g(x)] 



(ixx"[Ag(x) 



-l)"Ag(x)], 



(7) 
(8) 
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where q is the sum of the quarks with hehcity ahgned and anti-ahgned (g^ + q^) and Ag is 
the difference of the two ( q^ — q^)- On the lattice, this corresponds to the matrix elements 
of the operators 



111. ..fin 



(9) 
(10) 



respectively with D = ^{D — D) the difference between the covariant derivatives. Note 
that here the trace term which corresponds to the disconnected piece is not included in 
our calculation. Therefore, we will look at the difference between the up and down quark 
contribution where this contribution is negligible. 

In this work, we will only concentrate on the first moments of the unpolarized and helicity 
distributions, {x)u-d and {x)Au-d respectively. The corresponding renormalization constants 
are obtained using RI/MOM-scheme NPR as described in Sec. Ill A[ The renormalized 
moments are shown in Fi gur e |4] at each pion mass point. Again, we use help from chiral 



effective theory 



28 



29 



30 



3l| for these quantities to extrapolate to the physical pion point: 



{x)u-d 
{x)Au-Ad 



c 
c 



1 - 



2^1 + 1^2 



In 



mz 



(4vr/.)2 
(4vr/.)2' 



(11) 
(12) 



This chiral behavior is indicated in the blue line in Figure |H We see a strong curvature 
due to the chiral form; more light-pion points should be taken to eliminate extrapolation 
uncertainties. In the past, we have been finding these quantities to be about 50% larger 



than the experimental ones. (See Ref. ll|, [12] for example.) In this updated 2+ If DWF 
calculation, these moments are 0.133(13) and 0.203(23). Figure [5] shows a list of the latest 
calculations of the first moments of the unpolarized (left) and helicity (right) distributions. 
Here we can see the quenched or partially quenched approximation results either from our 

or from the QCDSF/LHPc|27, 



past DWF calculation pj, 



32 



seem to be consistent with the LHPC's mixed action calculation 



Our preliminary results 



23|, although more statistics 



in the near future will help us to clarify what role the staggered sea plays in these quantities. 
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FIG. 4: Renormalized first moment of the unpolarized (left column) and helicity (right column) 
distri butions, in terms of and their chiral ex trapolations 
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FIG. 5: Global comparison of the first moments of the unpolarized (left column) and helicity (right 
column) distributions, in terms of and their chiral extrapolations 

3. Transversity 

Another interesting quantity regarding the spin structure of nucleon is transversity. The 
moments of transversity are defined as 



5q 



dxx"[5g(x) - (-l)"5g(x)], 



(13) 



where 6q is the difference between the quarks with spin ahgned and anti-ahgned with the 
polarized target. On the lattice, this corresponds to the matrix elements of the operator 



'^/i//i.../^„ 



■n-l 



(14) 



Again, we only calculate the isovector quantity to eliminate the contribution from the dis- 
connected diagrams. 

We calculate the zeroth moment of transversity, (l)^^, and the results are given on the 
left-hand side of Figure [61 We observe rather weak dependence (linear extrapolation) on the 
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FIG. 6: Zeroth moment of transversity from our data: the chiral extrapolation to the physical pion 
point (left) and the global comparison among different lattice groups (right) 



quark mass. We use the chiral extrapolation formulation 



28 



29 



30 



3l| 



{x) 



5u—6d 



1 7^ — -^mt In 



2(47r^) 



+ e'(/i^ 



(4vrA) 



(15) 



and get 0.56(4) at physical pion mass point. This extrapolated value is significant smaller 
than the simulated pion mass point, which is of order 1 or so. We urgently need data from 
our lightest pion mass to confirm this rapid decreasing behavior. However, this is close to 



what has been found by LHPC with mixed action 



about 0.7; their data is listed on the 



right-hand side of Figure El Their results at each pion mass are consistent within our current 
statistics. Apparently, the total suppression of sea quarks plays an important roles, as seen 
in comparing our quenched and dynamical numbers; however, there is not much difference 
between 2 and 2+1 flavors. 



4- di 

The twist-3 first moment of the polarized structure function di is another interesting 
feature to consider. It is related to the polarized structure functions gi and g2 by the 
operator 

C?[3l] = ^'^l5[l3^4-lA'Ds]^p. (16) 

It mixes with the lower-dimensional operator O3I if the lattice fermions do not have chiral 
symmetry at finite lattice spacing; we are free of this problem in our calculation. Although 
it is not measurable in deep inelastic scattering of electrons on protons, it gives us some 
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FIG. 7: Bare twist-3 matrix element di 



expectation of the higher moment dn matrix elements. Figure [7] shows our isovector di 
matrix element results. We extrapolate the twist-3 matrix element to the physical pion 
mass and get d\^'^'^ = —0.002(2), which is consistent with zero. Combined with the small 
value of d2 from QCDSF33|], we conclude that the Wandzura-Wilczek relation between 
moments of gi and g2 3J], which asserts vanishing dn, is at least approximately true. 



III. HYPERONS 
A. Lattice Parameters 



We use a mixed action, meaning that the sea and valence fermions use different discretiza- 
tions. In our case, the sea fermions are 2+1 flavors of staggered fermions (in configuration 
ensembles generated by the MILC collaboration 35|]), and the valence fermions are domain- 
wall fermions. The pion mass ranges from 360 to 700 MeV in a lattice box of size 2.6 fm. 
The strange-strange Goldstone is fixed at 763(2) MeV, which unfortunately does not repro- 
duce the physical strange-strange Goldstone mass. The gauge fields are hypercubic-smeared 
to improve the chiral symmetry of the fermion, and the fermion field is Gaussian-smeared 
to improve the signal. The source-sink separation is fixed at 12 time units. Table [III lists 
details of the configurations in use. 



B. Strangeness Contribution 

Studying the strangeness content of the nucleon is important to understanding QCD. 
Since the nucleon has zero net strangeness, any contribution to nucleon structure observables 
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Label 


(MeV) 


niK (MeV) S conf. 


H conf. 


mOlO 


358(2) 


605(2) 


600 


600 


m020 


503(2) 


653(2) 


420 


436 


m030 


599(1) 


688(2) 


561 


561 


m040 


689(2) 


730(2) 


306 


319 



TABLE II: Configuration details 



is a purely sea-quark effect. Many experiments are devoted to understandin 



contributions to the elctromagnetic 
at JLab, SAMPLE at MIT-BATES 



brm factors of the nuc 



eon: HAPPEX 



40 



36 



strange quark 



37j and G0|38| 



41| . The experimental results 



3911, and A4 at Mainz 

reveal a small but non-zero strange contribution to the proton elctromagnetic form factors. 

To theoretically understand this nonperturbative physics, lattice QCD is a natural can- 
didate for applying first principles. However, to extract individual quark contributions, 
the difficultly of calculating disconnected diagrams (as shown in Figure [1]) must be taken 
into account. In the past, these diag rams have been directly calculated in the quenched 
approximation in lattice QCD HUH], giving values of ranging from —0.28(10) to 
+0.05(6). The Adelaide- JLab Collaboration used an indirect approach with the help of 
charge symmetry 4^ and chiral perturbation theory to correct for quenching effects, ob- 



taining —0.046(19) 431 . There does not seem to be consistency among these works. In this 
paper, we will use unquenched lattice data from mixed action with charge symmetry in the 
hope of bringing clarity to the chaos. 

From charge symmetry 42], 44 1 , one can derive the following relation of magnetic moments 
(/i^) of individual quark contributions (g^) and the disconnected contribution Ob for each 
octet baryon B: 



CdU^ + Cud'" + On] 
Crf^" + e^s" + Oh- 



The disconnected contribution to the proton Oat is 

1 1 - i?: 
3 



O 



N 



d 



M5 



(17) 



(18) 
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where i?^ is the ratio of the strange to down quark loop disconnected pieces, G\j/G'{j. 
Combining Eq. [T7] and Eq. [181 we get the strangeness magnetic moment from either 



G 



M 



.N 



2/iP + /i"-— (/i' 



/i 



(19) 



or 



M 



Rd 



u 



N 



2/i" ^(/i 



/i 



(20) 



The ^ (^) deviates shghtly from 1 due to SU{3) symmetry breaking. Using input from 
the well measured experimental quantities 45] fi^^ — fJ'^ = 3.618 or — fJ'^ = —0.599, 
we get a constraint on the stangeness content of the proton magnetic moments. The left of 
Figure [9] shows the range of G^^: it goes to negative values as indicated by the blue dashed 
line, while the positive ones lie on the red solid line on the ratio parameter plane. Note that 



since (2/1^ + /i") - (/i^"^ - /i^ 



0.057, one would need higher precision on the lattice data 



M- 



for ^ to make use of this constraint. Thus, in this work, we will only use Eq. |2D]for G 

We extrapolate to = using a dipole form for the magnetic form factor of octet 
baryons. Figure [8] shows our fuU-QCD lattice data on the ratios ^ (left) and ^ (right) 
as a function of m^. The physical limit is taken by naive linear extrapolation, since the 
ratio might cancel out higher-order dependence on pion mass in the chiral extrapolation; 
the two ratios are 1.03(13) and 1.04(13) respectively. We compare our extrapolated value 



43| in the left part of Figure O We 



with the calculation from Adelaide-JLab Collaboration 
find our statistical errorbar in S to be much larger than theirs. This might be because 
they use quenched configurations where the ensembles for different pion mass points are 
correlated; we extrapolated through different uncorrelated dynamical ensembles. Combined 
with Eq. [20l this gives a constraint (shown as a pink band in the right part of Figure [9]) 
on the GIj as a function of ratio of strange to up/down contribution to the proton i?^. 
Such a ratio will give smaller statistical error than calculating individual components, since 
fluctuations will be canceled in the ratio. However, at the moment, we do not have such 
a ratio calculated directly from the lattice. We quote the estimation i?^ = 0.139(42) |43| 
from chiral perturbation theory; this gives G%,j = — 0.066(12)stat(23)Rs with the dominate 
errorbar from the conservative estimation of R^, which is consistent with Adelaide-JLab 
Collaboration's number, —0.046(19). 
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FIG. 8: Magnetic moment ratios of proton to S and to H as functions of the pion mass 
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FIG. 9: (Left) Our to and u~ ratio (purple star) compared with Adelaide-JLab Collabora- 
tion's (green circle). The blue dashed line indicates negative values of G\,j, while the red solid line 
indicates positive ones. 

(Right) The pink band is the constraint for the proton strangeness magnetic moment from our 



data and the grey band indicates the given from Ref 



Similarly, with the help of charge symmetry, one can estimate the of the proton via 
the root-mean-square radius from the up quark connected diagram 46|: 



(21) 



The left-hand side of Figure [TO] shows the (r^)" at each pion mass and their chiral extrapola- 
tion according to Ref. 47|]. Note that this is the result from analyzing the number of config- 



uration listed in Table HTl only. Taking = 0.16(4) from chiral perturbation estimates 461] , 



we find G^^{Q'^ = 0.1 GeV) = 0.022(61), which is consistent with experimental values, as 
shown on the right of Figure [TOl The statistics will be further improved with the latest 
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region 



48t | with our preliminary result from mixed action and Ref. 



46| 



LHPC calculation. 



C. Hyperon Axial Coupling Constants 



The axial coupling constants, and (^se? have important applications such as in hyperon 



scattering and non- 



eptonic decays. Previously, there only existed predictions from chiral 



perturbation theory |49l] and large- A^c calculations: 



0.18 < -^ss < 0.36 
0.30 < < 0.55. 



(22) 
(23) 



Figure [TT] shows our lattice data and chiral extrapolation. Here we take a naive linear ex- 
trapolation against {rrij^/ f^^y. We find numbers consistent with the models, (7ss = 0.441(14) 
and = —0.277(11), but with much smaller errors. 



IV. CONCLUSION AND OUTLOOK 



In summary, we report the latest nucleon calculations using RBC/UKQCD 2+1-fiavor 
3 fm DWF ensembles with pion masses as light as 310 MeV. Our calculations show good 
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FIG. 11: Lattice data for <^SE and (7== and chiral extrapolation 



consistency with experimental values. Even those quantities, such as the first moments of 
the unpolarized quark distribution and helicity distribution, are chirally extrapolated to the 
physical pion mass in consistency with experimental values. We predict the zeroth moment 
of transversity, and we find the twist-3 matrix element d\ to be consistent with zero. 

In analyzing hyperon form factors, we use an indirect approach to get the strangeness 
of the proton magnetic moment, using mixed action. We find our dynamical result to be 
consistent with Adelaide- JLab Collaboration's quenched result (which used a chiral correc- 
tion for sea quark effects) and current experiments. The axial charge coupling of S and S 
baryons are also predicted with significantly smaller errorbars than estimated in the past. 
We will continue to increase statistics, especially in the light pion mass region to get even 
more accurate results in the future. 
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